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Abstract 

Brownian motion of particle interacting with atoms of ideal gas is discussed as a key problem of kinetics lying at 
the border between "dead" systems like the Lorentz gas or formal constructs of conceptual Boltzmannian kinetics 
and actual "alive" systems like mere gas possessing scaleless (1/f) fluctuations in their kinetic characteristics (e.g. in 
diffusuvity and mobility of the "Brownian particle"). 
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1. Introduction (once again about derivation of kinetics from dynamics) 

What is simpler than ideal gas? At least, when simply gas is too hard nut to crack? 

In this paper we want to consider Brownian motion of a particle interacting with infinite gas of atoms 
which do not interact between themselves. Our starting-point will be corresponding Bogolyubov-Born-Grcen- 
Kirkwood-Yvon (BBGKY) hierarchy of equations [1,2,3] for (n + l)-particle distribution functions of the 
"Brownian particle" and n atoms ( n = 0, 1, 2, ... ) or equivalent Bogolyubov functional evolution equation 
(BFE) for their generating functional [1]. Our task here is not to solve these equations but only discern them 
and discuss their possible solution in terms of reasonably introduced (n + l)-particle correlation functions, 
with taking in mind general invariance properties of BFE found in [4] and expressed by so-called "virial 
expansion of Brownian path probability distribution" and other "virial relations" . The latter were foreseen 
earlier in [5,6] and partly deduced in [7,8,9,10] from the generalized fluctuation-dissipation relations [11,12]. 

The modesty of our task is opposed to pretence of the legend existing among "nutcrackers" and stating 
that dilute gas (under the "Boltzmann-Grad limit" or, in other words, the "low density limit") is so much 
simple object that in respect to it the BBGKY hierarchy reduces to so-called Boltzmann hierarchy and 
then to the single Boltzmann equation or some its derivative like the Boltzmann-Lorentz equation [13]. 
The legend refers to works by Lanford on gas of hard spheres, e.g. [14,15] (see also [13] and references in 
[3] and [13,15]). However, careful reading of work [15] shows that factually it does not suggest a proof of 
the proclaimed statement. Moreover, the author himself indicates that already short time evolution (five 
times shorter than mean free-path time) reveals definite surfaces in n-particle phase spaces ( n > 1 ) where 
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n-particle distribution functions (DF) "do not converge to the desired products" of one-particle DF, thus 
preventing the desired proof. 

The mentioned (hyper-)surfaces are Vi — Tj = (vj — Vj)0 with r and v being particle coordinates 
and velocities and 9 arbitrary time parameter. That are surfaces made of various pre- or post-collision 
trajectories of gas particles. Hence, a correct conclusion what suggested itself in [15] was that generally 
statistics of collisions is not determined by one-particle DF, and BBGKY equations can not be reduced to 
Boltzmann equation (BE), even under the "low density limit". 

From physical point of view, this is quite trivial conclusion. Indeed, the "low density limit" in no way 
removes fluctuations in density of collisions (number of pairs in pre- or post-collision states per unit volume) 
as well as it does not remove fluctuations in density of particles (number of particle per unit volume) , sooner 
it strengthens all them. Therefore, even if local density of collisions was primitively thought as square 
function of local density of particles, an unbiased reasoning would result in the stated conclusion merely 
because of general inequality (A * B) 7^ (A) * (B) (with angle brackets denoting ensemble average) ! 

In later attempts to derive Boltzmann hierarchy from the BBGKY one, e.g. in [13], it was postulated that 
values of DF for an incoming, or pre-collision, state (at 9 < with above introduced 9 ) and responding 
outgoing, or post-collision, state (at > ) are equal. It is good idea (although wrong in literal quantitative 
sense). But elementary logics requires to apply it at once in all places of the BBGKY hierarchy where 
differential n-particlc Liovillc operators act on DF just realizing collisional transitions from in- to out- 
statcs. Instead, the authors applied their postulate in integral terms only, while in the differential terms 
neglected collisions at all! Of course, results of such arbitrariness hardly can be attributed to physics. 

More likely approach to kinetic theory of dilute gas (formally, in the Boltzmann-Grad limit) was suggested 
in [16] (or see [17]). There it was emphasized that reformulation of BBGKY hierarchy in terms of collisions of 
particles and "collision integrals" , in place of continuous interactions, is possible only if various consecutive 
stages of any collision process are presented in statistical ensemble under consideration with equal weights 
(probability densities). In other words, derivative of DF in respect to the "inner time of collision" 9 turns 
to zero inside a spatial region assigned to collision. This property never realizes exactly, but it must be 
postulated in order to construct a correct "collisional approximation" to BBGKY hierarchy. It means [16] 
that in any particular Liouville operator the replacement 



should be made within collision of i-th and j-th particles, where $(r) and — $'(r) = — <9$(r)/<9r are 
interaction potential and interaction force, respectively, and the first equality gives formal definition of the 
"inner time of collision" . To be more precise, the words "within collision" mean 



with To being characteristic interaction radius and A mean free path. 

Evidently, this ansatz by its nature is independent on shape of $(r) , therefore extends to hard spheres 
too. It is just what was deficient in speculations of [13] 1 . 

As the consequence, density of collisions determined by pair DF within collision drifts in space 
with velocity (vj + Vj)/2 . Since relative motion of colliding particles is included to collision, the latter as a 
whole moves with this centre of mass velocity! This is sufficient to forbid factorization of density of collisions 
into any quadratic functional (e.g. product) of two one-particle DF drifting with different velocities, Vj and 
Vj . Hence, pair DF of colliding particles is independent statistical characteristics of gas evolution. 



1 This loss demonstrates that hard spheres is treacherous model which may play a bad trick on its makers, at least when 
treated in the language of DF. One can see also what is principal error of [13,14,15] and similar mathematical "fantasies": 
although Lebesgue measure of regions "within collision" tends to zero under the Boltzmann-Grad limit, this is not a ground 
to cast such regions since just they determine actual evolution of the system! 
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Eventually, "collisional approximation" to kinetics of dilute gas results in an infinite hierarchy of kinetic 
equations [16,17,18] for usual one-particle DF and infinitely many less usual mutually independent n-particle 
DF. At n > 2 they represent ensemble-average densities of connected clusters of n — 1 (real or virtual) 
collisions. This new hierarchy reduces to Boltzmann hierarchy and then to BE in case of strictly spatially 
uniform statistical ensemble only (in this case, possibly (A * B) = (A) * (B) since ensemble average can be 
identified with infinite-volume average) . 

The word "kinetic" here means that interactions are represented by a ready-made "collision integrals" 
instead of potential $(r) . But, strictly speaking, this only gently simplifies the theory, since full list of argu- 
ments and parameters of n-particle DF ( n > 1 ) decreases by one only (namely, by the 6 ) [17]. Therefore, 
to derive benefit from the new hierarchy, some its further reasonable posterization and/or approximation is 
necessary. Two different approaches to this task were considered in [16,17] and [18,19]. 

Most principal result of [16,17] was presence of long-living statistical self-correlations in random ("Brown- 
ian") motion of any gas particle and corresponding low-frequency scaleless 1/f-type fluctuations ("1/f-noise") 
in its self-diffusivity (and mobility). In [18,19] this result was confirmed and seriously extended by showing 
that probability density distribution of diffusive path of a test particle possesses power-law long tails, natu- 
rally cut off at distances of ballistic flight (see also [9]). All this qualitatively justifies early phenomenological 
theory of 1/f-noise under realistic molecular Brownian motion [20,21] and, in turn, recently was justified in 
[7,8,9,10] and then in [4] basing on exact invariance properties of Lioville operator and BBGKY hierarchy. 

And all this is minimum of what is neglected by the Boltzmann hierarchy and Boltzmannian kinetics at 
all. Unfortunately, its modern admirers think on pioneer level with Boltzmann, as if they are unacquainted 
with Krylov's warning [22] that physical independence of particles or events on concrete phase trajectory 
(in real life) does not imply statistical independence of their images at statistical ensemble of trajectories 
("ensemble of lives"). May be, from here the secret wish arises to equalize in rights "dilute gas" and "Lorentz 
gas" (non- interacting particles in a lattice of fixed elastic scatterers). 

In fact, results of [16,17,18,19], as well as [7,8,9,10,4], say that "dilute gas" differs from "Lorentz gas" 
like a living matter differs from dead one. In the first, contrary to the second, random walk of a test 
particle never can be divided into statistically independent constituent parts, regardless of their durations 
and total observation time. In other words, from the point of view of time averaging, every trajectory of this 
random walk remains unique at arbitrary growing spatial-temporal scales, while in "Lorentz gas" almost all 
trajectories become identical (use of terms "ergodic" or "non-ergodic" I leave to mathematicians). A simple 
heuristic interpretation of this property was suggested in [18] and continued in preambles of [7,8]. 

Additional explanations can be found in [16,17,21] and in preambles or resumes of [9,19,20,23,24,25]. 

By the way, notice that the mentioned difference long ago is known as experimental fact concerning 
charge carriers in semiconductor crystals: the greater is concentration of hard immovable impurity atoms, 
the stronger is damping of relative 1/f fluctuations in carrier's diffusivity (and mobility) due to phonons (see 
e.g. [21] and references therein). 

The aforesaid makes clear urgency of Brownian motion of a particle interacting with ideal gas. Evidently, 
this problem mediates between "dilute gas" and "Lorentz gas" , especially under the Boltzmann-Grad limit 
when radius of the Brownian particle is arbitrary small as compared with its (fixed) mean free path. We 
expect that Brownian motion in ideal gas also is alive, contrary to dead random motion among static 
scatterers. Let us consider this expectation starting from [4] and taking in mind experience of [16,18]. 

2. Functions and equations of the model 

As in [4], in the meantime we confine ourselves by Browniam motion in equilibrium fluid (to be ideal gas 
here), for concreteness, in the framework of canonical Gibbs ensemble of initial conditions for fluid atoms. 
What is for the Brownian particle (BP), again initially it is certainly located at given point Ro . But its 
initial momentum is random and obeys equilibrium Maxwell distribution (with the same temperature T ) . 

Let $(r) denotes potential of interaction between BP and atoms. Since now atoms do not interact 
one with another, the relations between DF F n (t, R, r("), P, |R ; v) and correlation functions (CF) 
V n {t, R, P,p(™)|Ro; v) , with = {ri... r„} , p(") = {pi... p„} and v standing for mean density of 
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gas (density at infinity), the BBGKY equations and the Bogolyubov functional evolution equation (BFE) 
all strongly simplify. The BBGKY equations take form 



dt M <9R dP ' ' 



dF d { 

—^ = L n F n + V — J $'(R-r n+1 )F n+1 

ra+l 

at n > 1 , with J k ... = J ... dp k dr k , 
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dP dp. 



(3) 



and obvious initial conditions: 



F n (t = 0, ... ) = 5(R - R ) G M (P) J] G ™(Pj) E{v 3 ) , 

G m (p) = (2^Tm)- 3 / 2 exp(-p 2 /2T™) , £(r) = exp[-$(r)/T] 



(4) 



The connection between generating functional of DF 

oo 

T{t, R, P, V |Ro; v} = F (t, R, P|R ; v) + V - 



n=l 



F n (i, R, rW, P, P (n) |R ; v) Y[ V(r j5 Pj ) 



and quite similarly defined generating functional of CF [4], V{t, R, P, ip |Rq; , becomes 



!F{t, R, P, ^ |Rq; v} = exp 



i/ J G m ( P )£(r-R)V(r,p)d P dr 



V{t,R,P, V|R-o; 4 



= exp 



j/ / £(r - R) 0(r) dr 



(5) 



V{t,R,P, V|Ro; f} , 

where </>(r) = J ^(r, p) G m (p) dp . Recall that in fact this is definition of CF. According to it, 
F (t,R,P|R ;^) = Vb(t,R,P|R ;^) , 

Fi(t,R,ri,P,pi|Ro;i/) = V r o(t,R,P|R ;i/)S(ri -R)G m (pi) + ^i(t,R,ri,P,pi|Ro;i/) , 
and so on. The BFE, that is compact functional form of BBGKY hierarchy, now reads 
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with .t = {r,p} and J ... — J ...dx — J ...dpdr, thus turning into first-order differential equation in 
respect to ip(x) . Substitution of (5) to (7) yields equivalent BFE in terms of CF: 
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with E'(r) = dE(r)/dr = -$'(r)£(r)/T . 



Variational differentiations of (8) return us to BBGKY equations rewritten in terms of CF: 
dV P dV Q 



dt M <9R 



"dPj *'( R - r i)^ . ( 9 ) 
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BV B f ( P B \ 

= LiVi + v-Qpj *'(K-r 2 )V 2 + rG ra ( Pl )£'(n-R)^- + — jVo , (10) 

2 
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(11) 



V n - 1 (...r j _ 1 ,r j+1 ... Pj _ 1 ,p j+1 ...) 



Initial conditions and boundary conditions to these equations and (8) are very simple: 

V n (t = 0,...) = S n , <5(R-R )Gm(P) , V{t = 0, R,P, V|Ro;^} = ^(R-Ro)Gm(P) (12) 
(I have to detect misprint in similar formula (13) in [4]: there factor Gm(P) was missed!), 

K>o( -rfe ...)-» at r fe -R^oo (13) 

At last, consider the invariance group found in [4] (a group of such transformations of arguments of 
generating functional determined by (8)-(13) which do not change value of the functional). In case of ideal 
gas it simplifies to 

V{t, R, P, a + V | Ro; y] = V{t, R, P, V/(l + a) | R ; (1 + a) v} , (14) 

where cr(r, p) = a = const is arbitrary constant from interval —1 < a < oo . This functional identity implies 
exact "virial expansions" of CF: 

K(i,R,r (n) ,P,p (n) |Ro; v + trv) = 

= K(i,R,r("),P,p(«)|R ; «/) + f] j ... j K+fe(i, R, r^ n+k ^ , P, p("+ fe ) |R ; „) ^ 

k = 1 n+l n+k 

Their infinitesimal form yields exact "virial relations" 

dK(t,R,r("),P,p(»}|R ;^) f (n+1) (n+1) 

n+l 

At n = , in particular, we obtain connection between density derivative of the probability distribution of 
path, AR = R R , of the Brownian particle and integrated pair CF: 

aVo(t,R,P|R ; v) /■„,.„ „ m v, , 

— = / ^i(i,R,r,P,p|R ; i/) drdp , (17) 

where, of course, both sides depend on current position of BP, R , and it starting position, Ro , through 
their difference AR = R — Ro only. 

Relations like (17), as combined with identity (6) and trivial positivity of DF Fi , lead to principal 
consequence [4,7,8,9,10] mentioned in Introduction: asymptotic of BP's path distribution, 

V (t,AK;u) = yV (t,R,P|R ; V ) dl> , (18) 

at t » r , with t being mean free-path time of BP, as a function of AR has power-law long tails lasting 
up to |AR| - ty/T/M . 



Such statement strikingly contradicts Boltzmannian kinetics which always gravitates towards the "law of 
large numbers" and Gaussian asymptotic characterized by short exponential tails. To avoid new repeating 
myself, I ask dear reader to search for proper comments in [8,9,10,16,17,18]. But one not popular truism 
deserves repetition: in statistical physics, unlike probability theory, one should not rely on the "law of large 
numbers" since practically independent events or quantities may be nevertheless statistically dependent. It 
would be a good thing to understand in detail how all this realize in case of Brownian motion in ideal gas. 

Notice that due to simplicity of BBGKY equations (l)-(2) or (9)-(ll), in comparison with general case 
[4,16,18], one can easy verify virial relations (16) by deriving them directly from (9)-(ll). 

3. Relative coordinates and characteristic function 

Eventually, we would like to obtain the probability distribution (18) of BP's path or its characteristic 
function, that is Fourier transform 

V (t,ik; v) = J exp (ik ■ AR) Vo(t, AR; v) dAJH (19) 

Therefore, first, let us consider all the DF and CF as functions of AR = R Ro and relative distances 
Pj = r j — R . Such change of spatial variables implies operator changes 

d d >— ^ Odd 

dR. ^ MR ~ ^ dp~ ' dr~ ^ ~dp~] 

everywhere in (l)-(2), (9)-(ll). Second, make Fourier transform in respect to AR and go to functions 

V n (Mk,p< n \P,p( n );i/) = y"cxp[ik- (R-Ro)] y„(t,R,R + p (,l) ,P,p (ll) |Ro; v) dH (20) 

Third, introduce velocities V = P/M and Vj = Pj/m of BP and gas atoms, and new operator 

d . / d d 



Then BBGKY equations (9)-(ll) take the form 
dV .„ „ x „ d 



0j i(k-V)V - u — J &{Pi)Vi , (22) 

l 

^ = i(k-V)Vi + A(l)Vx -v-HpJ &(P2)V 2 + G m ( Pl )E'( Pl ) (v + T-^V , (23) 

2 

= i(k-V)V n + ^>(j)K - u — J $'(p„+i)K+i + 



j- 1 n+l 

OP 



+ 5^G m (pj)^'(Pj)' (v + T^)v r n _i(...p j _i,pj + i...pj_i,pj+i...) 



(24) 



with initial and boundary conditions as follow, 

V n (t = 0,...) = 5„,oG M (P) , K>o(-Pfc-) -> at p fe -» oo , (25) 
and BP's path distribution (18) under interest presented by 

rfk 

(27p ' 

where d is space dimension ( d = 3 by default). 



Vb(t,AR; v ) = J cxp (-ik • AR) ^ V (Mk,P; i/) dP , (26) 
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Taking into account quite obvious identity (V + T d/dP) Gm(P) =0 it is easy to see that at k = 
the system (22)- (24) always stays in its initial state: 



V (t,ik = 0,P; v) = G M (P) , 



K(Mk = o,p (n) ,P,p (n) ; v) = /"v n (t, R,R + P (n) ,P,P (n) |Ro; v) dn = (n>0) 



(27) 



First of these expresses mere normalization of BP's path distribution and besides says that (unconditional) 
probability distribution of BP's velocity in equilibrium gas always stays equilibrium. 
Collecting all CF V n (t, ik, p( n \ P, p(™); v) into generating functional 



V{t, ik, P, V ; v} = V (t, ik, P; v) + J2 



n=l 



V 



V„(t,ik,pW,P,pW;i/) JJ <p(Pj,Pj) 

3=1 



one may replace all equations (22)-(24) by single functional equation: 
dV 

m= £v 

with evolution operator 



£ = i(k • V) 



ijj(x) 



d_ 

d~ P 



(-- 


-V 
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9p y 


Sip(x) 



(28) 



(29) 



G m (p)E'(p)iP(x) 



9P 



and initial condition 

V{t = 0,zk,P,^;4 = G M (P) 



(30) 



In essence, of course, this is full equivalent of (8) following from (8) under change tp(r,p) — ► ^>(R + p, p) 
and correspondingly 



^ 9R 



dp / dp 



dip(p,p) 



dp 5ip(p,p) 

4. Boson representation and path integrals 

1. According to formulas of previous section, direct formal exact solution of BBGKY equations (22)-(24) 
in respect to the characteristic function (19) of BP's path distribution can be written as 



V Q (t,ik-v)= jv{t,ik,P,i> = 0;v} dP = J dP e t£ G M (P) 



V> = o 



(31) 



where the evolution operator £ represents, in respect to BP's momentum, a linear combination of two 
vector operators, V + T9/9P and —Td/dP (factor V in first row of (29) is their sum). Let us normalize 
them so that their components commute one with another exactly as boson birth and annihilation operators: 



At 



Vtm 



d_ 
dP 



A = 



V + T 



d_ 
dP 



(32) 



Then notice that AG M (P) = and J dP At ... = . Therefore the Maxwell distribution, G M (P) = 0, 
can be treated as ket-vector of ground state of "quantum harmonic oscillators" represented by A and At , 
while integration over momentum, J dP ... , acts as multiplication by bra- vector of the ground state. 
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Besides, notice or recall that operators of multiplication by ip{x) and differentiation 5/5^p{x) also behave 
as boson birth and annihilation operators, and we can write 



a)(x) = c(x)ip(x) , a(x) = c 1 (x) 



8ip{x) 



a(x)a\y) - a* (y)a(x) = 5(x - y) 



(33) 



with arbitrary fixed c{x) ^ . At that, unit from the right of them and putting on tp — after their action 
again represent ket- and bra-vectors of vacuum state, respectively. A proper choice of c(x) is 



c(x) = ^ V G m (p)E{p) 

For further let us introduce also vo = \JTJm , u a = ^/TjM . 
Then after some algebra formulas (29) and (31) can be rewritten as 



V (t,ik; v) 



, t£ 



10) , 



£ — £\ + £2 + £ 



(34) 



3 , 



where |0) is common ground state ("vacuum state"), and three parts of the evolution operator are linear, 
quadratic and cubic forms of the birth and annihilation operators, respectively: 



£ l = i(k- (A + At))« 
d 



£2 = 



a\x) 



a(x) + u / c(x) 



■ [a(x)A^ - a\x)A] , 



= u Q J a){x) 



a(x) 



(35) 
(36) 

(37) 



In such way calculation of exact characteristic function of BP's path transforms to calculation of the 



"vacuum- vacuum transition" amplitude, 



- an 



|0) , for excited system of interacting bosons, or quan- 



tum oscillators, with cubic Hamiltonian 7i = i£ . Due to the boundary conditions (13), (25) one can think 
that <v(x)a(x) — > at p — > 00 in those sense that far enough oscillators almost surely stay in their ground 
states. Therefore H can be treated as Hermitian operator. 

2. On this way, next possible step is use of so-called holomorphic form of path integrals (see e.g. [26,27,28]). 
According to it, we can replace (34) by path integral 



V (Mk; v) = (0| e C£ |0) - J exp J J 



I {A* A - A* A) + £(A*,A) 




dA*dA 
2iri 



(38) 



supplemented with edge conditions A*(£ = t) = , A(£ = 0) = , where A* = A*(0 = {a*{x, £), A*(^)} 
and A — A(£,) — {a(x, £), A(£)} are holomorphic images of all above introduced birth and annihilation 
operators in the form of complex- valued time functions, and the dot means derivative in respect to time £ . 
Integrating (38) first over variables a* (#,£), a(x,£) deputed by aJ(x) and a(x) one obtains 



V (t,ik; v) = J exp <| J 



(A* • A - A* • A) + t(k • (A + A*)) u 



x exp 



J A*(&) G^{6,6, A*, A} Ap{(a) d^d^ 



*>a>6>o 



n 



dA*dA 

{2iTi) d 



(39) 



with edge conditions A*(t) = , A(0) = , where repeated indices imply summation, the kernel G ai/ g is 
functional of A*(£) = and A(^) = expressed by 
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GapiSuh, A*, A} = || J c{x)& a { P ) exp | A(A*(0,A(0)d^ 



and A(A*, A) is sum of evolution operators from (36) and (37), 



(40) 



A(A*,A) = -v~ + + 



(A* + A).- + ( A*-A).^ 



(41) 



(clearly, this is image of the Liouville operator (21)). Here and below exp designates chronologically ordered 
exponential and, as before, / ... = J / ...dp dp. With use of obvious identities 



A(A*,A)c(x) = - 



mo $'(p) ■ Ac(x) 



, A T (A*,A)c(x) = 



u o ■ A* c(x) 



(where T symbolizes transposition) one can transform (39) into 

V (t,ik;v) = yexpN 1 (A* • A - A* • A) + i(k • (A + A*)) u 



x exp < / c(x) exp 



A(A*,A)d£ 



Lo 



i n 



X 

dA*dA 

(2ttz) 



(42) 



Exact analogue of formula (42) was obtained in [29], in slightly different notation, by means of so-called 
"stochastic representation of deterministic interactions" (see also references from [29]). 
Integration of (38) at first over A and A* results in another kind of holomorphic path integrals: 



Vb(t,ik; v) 



exp 

ii 



— (d*a — a* a ) + a* Aq a 



dxd£> x 



x exp jug f d£i f dfc [ik + K*(a*(ei),o(^i))]-[*k + K(a*(6),o(6))] I JJ" 
loo J 

where edge conditions a*(x,£ = t) = , a(x,£ = 0) = must be satisfied, 

A = -v.- + $' (p) .- 

is Liouville operator of atom interacting with BP (fixed at coordinate origin) , and 

d <P'(p)~ 



da* da 
2iri 



K*{a*,a) = - J c(x)-^ a*(x)dx + J a*(x) 



K(o*,a) = / c(x) a(x)dx + / a* (a;) 
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a(x) (ix , 



dp 2T 



a(x) dx , 



(43) 



where j ...dx — j j ...dp dp . 

Unfortunately, a correct non-perturbation treatment of such strongly non-Gaussian path integrals as (39) 
or (42) or even (43) is in itself non-trivial problem [27]. To make use of them, one has to construct some 
perturbation expansion (e.g. over m/M in the limit of hard BP). 



5. Generalized boson representation and continued fractions 



A different boson-like representation is prompted by characteristic triple-diagonal form of the system of 
equations (22)-(24). Considering V n as components of vector in the Fock space, let us define in this space 
birth and annihilation operators, ft and A, a* and a, as follows: 



P + TM 



dP 



G m (p„+i 
■ a t V n , 



dP 



AV n>0 = -v— ^'{ Pn )V n = —-aiV n , AV = 0, slV = 



<9P 



where symmetry of all the CF V n is taken in mind. According to this definition, 

{a a a) } - a^a a ) = a 2 5 a p , 



(A4 f - ftA)V n = a 2 



i =1 



(44) 



(p)dp , 



if we assume that is spherically symmetric potential and define projection operator II as 



U(j)V n = dG m ( Pj )E'( Pj )- / &{p 3 )V n d Pj dv 



J ^{P)-E'(p) 



dp 



Besides, introduce operator A k by equalities A k Vo = iik ■ V) V a , 



A k y„ 



i(k-V) + ^A(i) 



Then formal solution of equations (22)-(24), under initial conditions (25), in respect to Laplace transform 
of Vo(t, ik; v) , can be represented by continued fraction: 



oc 

/ 



e' zt V (t,ik; v)dt = 



I 



dP 



dP 



z-Ai 



Ai 



z - Ai 



'ft 



ft 



G M (P) 



ik - V 



d 



— T^(z,ik){P^ + TM 



G M (P) 



(45) 



where T^p begins recursive chain of operators 
T™(z,ik) = a c 



(46) 
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6. Stochastic form of the boson representation 



Let A^ , A be a set of pairs of boson birth and annihilation operators, such that A a Ap — ApA a — 5 a p . 
Then any path integral like (38) or (39) or (42) or (43), 



J expi^j ±(A*A - A*A)dd + ¥{A* (£M(0> j II 
= J exp[W{A*{£),A(0} 1 dM{A*,A} , 



dA*dA 
2-ni 



with arbitrary functional exp [F{*4*(£), A(£)} ] denned on interval < £ < t, can be formally considered 
as averaging of this functional over Gaussian "probabilistic" measure dM{A* ,A} : 

J exp [ ¥{A*(£),A(Z)} ] dM{A*,A} = (exp [ F{yl* (0,-4(01 ] ) (47) 
The measure dAA{A*,A} is completely characterized by corresponding pair correlation functions: 

(AcWApfa)) = (A* a (ti)A* (t 2 )) = , (A a (h)A}(t 2 )) = 5 a/J e(ti-t 2 ) , (48) 

where Q(t) is Heaviside step function. To prove these equalities, it is sufficient to calculate Gaussian integral 
which represents characteristic functional of the stochastic processes A* (t) and A(t) , 



(t \ r t t' 

J [b(0-A*(0 + b*(0-AZ)}(% \dM{A*,A} = exp J dt' J dt" b*(t')-b(t' 



, (49) 



where the mentioned edge conditions are taken in mind, and b(t) , b*(t) are arbitrary time functions. 
In this sense, in particular, 

V (t,ik; v) = (0| e t£ ^'*> |0) = ^cxp j £(A'(Q,A{0)d^ 

with operator function £{A\ A) defined by (35)-(37), and 



(50) 



V (t,ik; v) = ( exp 



* 1 \ 

i«ok-y"[A(0+A*(0]de- / / A*(a)-GUi,6,A*,A}-A(£ 2 ) \ 



(51) 



*>6>?2>o 

instead of (39), with Gaussian random processes A*(£) and A(t) defined by 

(A a (£ 1 )A /3 (i 2 )) = ( A* (ii) A^(i 2 ) ) = , ( A a (ti) A* (t 2 ) ) - 5a/»e(ti-t 2 ) 

Such "stochastic" point of sight at calculation of path integrals can make it more constructive. 
Notice that for any functional F = F{A*(£,) 7 A(£,)} arranged like second term in (51), 



(52) 



J ^(ajG^te-^.^.^^K^Wi, (53) 
*>6>«2>o 

where G a p{^i,^2, A*,A} involves -4*(£) and A(£) from interval £2 < £ < £i only, the identities 

(¥ n {A*(0,A(0}) =0 (n>0) , (exp[FM*(0M(0}]> = 1 
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take place, that is by itself such functional is identical to zero. In special case of quadratic (bilinear) func- 
tional, when <G a/3 {£i,£ 2 , A*, A} = G a p(£i - 6) , 



/ 



exp 



*>£l>?2>0 



t' 



(54) 



= cxpi J dt' J dt" b* a (t')Q af} (t' -t")b {t") 





where matrix function Q a p{t) is defined by 



oo oc 

q = e{i + ® g ® e}- 1 , J e - zt Q(t)dt = z + j 



i e - zt g(t)dt 





(55) 



with symbol <g> standing for time convolution. The first of the latter expressions remains valid also when 
Gap — Ga0(£i,&) is not a difference kernel. 

7. Discussion and conclusion 

At present, unfortunately, none of formally exact expressions (31), (34), (39), (42), (43), (45) or (51) 
can be calculated exactly or at least correctly, at least in the long-time limit t/r p oo with r p denoting 
relaxation time of BP's momentum. Therefore it remains only to discuss the convenient approximation of 
exact theory and establish discreditable invalidity of this approximation from the point of view of exact 
"virial relations" (15)-(17) (see [4,5,6,7,8,9] and remark at end of Sec. 2). 

1 . The mentioned approximation follows from the chain (22)- (24) if we cut off it already at second level, 
i.e. neglect second-order (three-particle) correlation and thus all higher-order correlations. This is just what 
one always makes (knowingly or unknowingly) when creating Boltzmannian kinetics. Then 



fi(Mk,P,p,p; v) = J e [ lk - v + X ]( t - t ')G m (p)i?'(p)- (v + T-^)v (t',ik,P- u)dt' 
o 

and the first BBGKY equation (22) turns into closed kinetic equation 



(56) 



V (t') dt' 



(57) 



with Vo(t) = Vo(t, ik, P; v) on the right and operator- valued kernel 

G a p(9,ik) ^^Jdpjdp $'Jp) e^ +1 ^ e ^(p)G m (p)E(p) 



(58) 



From the point of view of continued fraction (45), this is "one- loop approximation" when one substitutes 
zero for T^i and thus leaves two floors of the fraction only. 

Evidently, this kernel is a sharp function of 9 , with width nearly equal to time duration of BP-atom 
collision, To = ro/vo. At that, principally we are interested in the long-time long-range limit only, when 
t/r p — > oo and k — > under k 2 t — const . Therefore factor i(k ■ V) in (58) can be neglected, and equation 
(57) reduces to the "Boltzmann-Lorentz equation" 



^ = i { k.\)V Q + BV Q , A 



(59) 
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where B plays role of linearized collision operator. 

Interestingly, above "derivation" of Boltzmannian kinetics had not required the Boltzmann's "Stosszahl- 
ansatz" . Though the latter is necessary if one wants to transform B into standard "collision integral" . This 
observation shows that the heart of Boltzmannian kinetics is neglect of third- and higher-order correlations 
(and thus, in essence, neglect of true second-order correlation). 

Recall that A contains P and d/dP , therefore in general operator G(9,0) significantly depends on P 
and d/dP . But under the long-range limit it effectively interchanges to J dP Q(6, 0) Gm(P) ■ Therefore, 
regardless of details of B , long-range asymptotic what follows from (57) or (59) is the Gaussian one: 

V (t,ik,P;i/) -> cxp [-£>(!/) k 2 t]G M (P) (l + Tf ik-P + ...) , , ^ 

v M / (60) 

V (t,i\s.;v) -> exp[-D(v)k 2 t] , 
where BP's diffusivity D(y) and momentum relaxation time are presented by 
D(v) = u t p cx v~ , 

oc 

Jdejdpjdpjdp • e 2e $'(p) G m (p)G M (P)E(p) = (61) 

M + m V J /' V_V| S d v - V D G ™(P) G M(P)dpdP 

In (60) last multiplier of V a (t,ik,P;v) is important for (56), while the dots replace unimportant terms, 
and in (61) £ is full effective cross-section of BP-atom collisions. 

All that seems beautiful till one confronts equations (56) and (59) to the simplest virial relation (17). The 
latter requires that 

dV (t,ik; v) 



v 

TMd 

o 

2m 



dv 



J J yVi(t,ik,P,p,p;i/)dpdpdP (62) 

Combining this identity from exact theory with (60) and (61) we see that if the mentioned approximation 
was correct then we would have 

JJJ Vi(t,ik,P,/>,p; v) dpdpdP -> D ^ kH e -D(,)^t (63) 

In fact, however, expression (56), as combined with (60) and (61), after quite standard (although rather 
troublesome) manipulations yields 

/ / / Vl (*' ik ' P '^'P; ")dpdpdP^ 

^ _^j exp - l - u ltf(t-t'f-D{v)k 2 t' (t-t') dt' -> -i e - D H k2 * , ^ 



since exponential in (56) corresponds to free ballistic flight of BP after its single collision with an atom. 

The difference between (63) and (64) is qualitative, and it says that the conventional picture, including 
the Boltzmann-Lorentz equation, is far from truth! 

Such strong discrepancy came from our neglect of the third-order correlations in (23), i.e. V^'s contribution 
to (23), and thus neglect of all higher-order correlations. However, the "virial expansions" (15), in particular, 
first of them, as written via Fourier transforms (20), 



V (Mk,P; v + av) = V (t,ik,P; v) + J ... J V n (t, ik, p™, P, p«; v) 

71 = 1 1 n 

(recall that a > — 1 ) , help us to understand that cutting of any even high correlations is bad idea. 



(65) 
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Indeed, gas densities on two sides of (65), v + uv and v , can be different in arbitrary strong proportion 
< a + 1 < oo . Therefore any cutoff in the infinite series in (65) would give a faulty result (like e.g. cutoff 
in the series representing exponential function). Hence, in practice all terms on right-hand side of (65) are 
equally important, regardless of value of gas density! 

Consequently, all equations of the hierarchy (9)-(ll) or (22)-(24) are equally important for correct analysis 
of the BP's path probability distribution, even in the "low density limit" ("Boltzmann-Grad limit"). This 
means that Boltzmannian kinetics is not a true "zero-order approximation" of rigorous theory in respect to 
the density. Or, better to say, true kinetics has no literal "zero-order approximation" at all. 

2 . The reason for all this was explained more than once in [20,21] twenty five years ago, then in [16] 
and later in [17,18,19,23,24,25] and [6,7,8,9,10] (and in principal sense anticipated in [22] sixty years ago). 
Indifference of many-particle system to a number of happened events of definite kind (BP's collisions with 
gas atoms here) leads to scaleless fluctuations ("1/f -noise") in "mean number of events per unit time" and 
related quantities (e.g. BP's diffusivity and mobility here). 

In spite of this understanding, one can envy creators of conventional kinetics (see Introduction): it already 
resolved all its problems. But, from the other hand, this is kinetics of an invented tiresome and "dead" 
world. Our consideration demonstrated that even such simple world as Brownian particle interacting with 
ideal gas apparently is "alive" and interesting. Nobody is able to predict what "number of collisions per unit 
time" will meet this particle in particular life. Hence, real theory is not in the past, it is yet in the future. 
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